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Abstract . This report is the second of a series dealing with _£he 
determination of optimal trajectories for the aeroassisted flight experiment 
(AFE). The AFE refers to the study of the free flight of an autonomous 
spacecraft, shuttle-launched and shuttle-recovered. Its purpose is to gather 
atmospheric entry environmental data for use in designing aeroassisted 
orbital transfer vehicles (AOTV). 

It is assumed that: the spacecraft is a particle of constant mass; 
the Earth is rotating with constant angular velocity; the Earth is an oblate 
planet, and the gravitational potential depends on both the radial distance 
and the latitude; however, harmonics of order higher than four are ignored; 
the atmosphere is at rest with respect to the Earth. 

Under the above assumptions, the equations of motion for hypervelocity 
atmospheric flight (which can be used not only for AFE problems, but also 
for AOT problems and space shuttle problems) are derived in an inertial system. 
Transformation relations are supplied which allow one to pass from quantities 
computed in an inertial system to quantities computed in an Earth-fixed system 
and viceversa. j 

Keywords . Flight mechanics, hypervelocity flight, atmospheric flight, 
coordinate systems, equations of motion, transformation techniques, optimal 
trajectories, aeroassisted flight experiment, aeroassisted orbital transfer, 
space shuttle reentry. 
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1 . Introduction 

This report is the second of a series dealing with the determination 
of optimal trajectories for the aeroassisted flight experiment (AFE). The 
AFE refers to the study of the free flight of an autonomous spacecraft, 
shuttle-launched and shuttle-recovered. Its purpose is to gather atmospheric 
entry environmental data for use in designing aeroassisted orbital transfer 
vehicles (AOTV). 

It is assumed that: (a) the spacecraft is a particle of constant mass; 

(b) the Earth is rotating with constant angular velocity; (c) the atmosphere 
is at rest with respect to the Earth; (d) the Earth is an oblate planet, 
and the gravitational potential depends on both the radial distance and the 
latitude; however, harmonics of order higher than four are ignored. 

Under the above assumptions, the equations of motion for hypervelocity 
atmospheric flight (which can be used not only for AFE problems, but also 
for AOT problems and space shuttle problems) are derived in an inertial system. 
Transformation relations are supplied which allow one to pass from quantities 
computed in an inertial system to quantities computed in an Earth-fixed system 
(Ref. 1), and viceversa. 

Previous Research . Previous research on the topics covered here can 
be found in Refs. 2-12. For the general theory of flight paths and coordinate 
systems, see Refs* 2-3; for the equations of flight over a spherical Earth, 
see Refs. 2-4; for the perturbed motion about an oblate Earth, see Ref. 5; 
for AFE problems, see Ref. 6; for reentry problems, see Ref. 7; for methods 
of orbit determination, s,ee Refs. 8-9; for the values of the astrophysical 
quantities* see Ref. 10; for the values of the characteristic constants of 
the oblate Earth, see Refs. 11-12. 
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Outline . Section 2 contains the notations, and Section 3 defines the 
basic coordinate systems. The relations between coordinate systems are 
discussed in Section 4, and the angular velocity (or evolutory velocity) 
is introduced in Section 5. The kinematical equations for an inertial system 
are derived in Section 6, and the dynamical equations are obtained in 
Section 7. Section 8 summarizes the results, and Section 9 presents the 
transformation relations which allow one to pass from quantities computed in 
an inertial system to quantities computed in an Earth-fixed system, and 


viceversa. 
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2. Notations 

Throughout the paper, the following notations are employed: 
a = acceleration, m/sec^; 

A = aerodynamic force, N; 

C D = drag coefficient; 

C L = lift coefficient; 

Cg = side force coefficient; 

D = drag force, N; 

f = latitudinal component of the gravitational acceleration, 

2 

m/sec ; 

g = radial component of the gravitational acceleration, 

2 

m/sec ; 

J 2 = characteristic constant of the Earth's gravitational field; 

Oq = characteristic constant of the Earth's gravitational field; 

J 4 = characteristic constant of the Earth's gravitational field; 

L = lift force, N; 

m = mass, kg; 

M = Mach number; 

Q = side force, N; 

r = radial distance, m; 

r g = equatorial radius, m; 

r = polar radius, m; 

P H 

R = Reynolds number; 

2 

S = reference surface area, m ; 

T - thrust force, N; 

2 2 

U = Earth's gravitational potential, m /sec ; 
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V = velocity, m/sec; 

W = gravitational force, N; 
x = Cartesian coordinate, m; 

y = Cartesian coordinate, nr, 

z = Cartesian coordinate, m; 

a = angle of attack, rad; 

Y = path inclination, rad; 

6 = longitude, rad; 

y = bank angle, rad; 

3 2 

y = Earth's gravitational constant, m /sec ; 

3 

p = air density, kg/m ; 

0 = sideslip angle, rad; 

4> = latitude, rad; 

X = heading angle, rad; 

u = angular velocity of the Earth with respect to an 
inertial system, rad/sec; 

u . .■ angular velocity of the local horizon system with respect 
hi 

to the inertial system, rad/sec. 

Subscripts 

b = body axes system; 
e = Earth axes system; 
h = local horizon system; 

1 = inertial axes. system; 

w = wind axes system. 

Superscripts 

• = derivative with respect to time; 

-* = vector quantity. 
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3. Basic Coordinate Systems 

The basic coordinate systems for flight over a spherical Earth are the 
inertial axes system Ox.y^., the Earth axes system Ox e y e z g , the local 
horizon system Px h y h z h , the wind axes system Px^z^ and the body axes 
system Px b y b z b . 

3.1. Inertial Axes System . The inertial axes system Ox^z^ is a 
Cartesian reference frame defined as follows. Its origin 0 is the center of 
the Earth; the z^-axis is aligned with the axis of rotation of the Earth 
and is positive northward; the axes x.»y^ are orthogonal to the z.-axis 
and are directed radially; the trihedral Ox.y.z^ is right-handed. In 
particular, the plane x^ ,y^ contains the fundamental parallel (the Equator); 
and the plane x.,z. contains the fundamental meridian (the Greenwich meridian) 
at a particular time instant, the time instant t = 0. The symbols 

denote the unit vectors of the inertial axes system; these unit vectors are 
time invariant by definition. 

3.2. Earth Axes System . The Earth axes system Ox e y 0 z e is a Cartesian 
reference frame which is rigidly attached to the Earth. Its origin 0 is 
the center of the Earth; the z e ~a*is is aligned with the axis of rotation of 
the Earth and is positive northward; the axes x g ,y e are orthogonal to the 

z -axis and are directed radially; the trihedral Ox y z is right-handed. 

0 C C C 

In particular, the plane x e ,y g contains the fundamental parallel (the Equator); 
and the plane x g ,z e contains the fundamental meridian (the Greenwich meridian) 
at all time instants. The symbols i' e »J e »k e denote the unit vectors of the 
Earth axes system; these unit vectors are time dependent. 

3.3. Local Horizon System . The local horizon system Px h y h z h is a 
Cartesian reference frame defined as follows. Its origin P is identical with 
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the instantaneous position of the spacecraft; the z^-axis is directed 

radially (that is, vertical) and is positive downward; the axes x h ,y h 

are orthogonal to the z h -axis (therefore, they are tangent to the spherical 

surface through P; they form the so-called local horizon plane); the 

trihedral Px h y h z h is right-handed. In particular, the x h ~axis is tangent to 

the local parallel through P and is positive eastward; the y h ~axis is tangent 

to the local meridian through P and is positive southward. The symbols 

i i_ * j ^ * k. denote the unit vectors of the local horizon system, 
h h h 

3.4. Wind Axes System . The wind axes system P*^^ is a Cartesian 

reference frame defined as follows. Its origin P is identical with the 

instantaneous position of the spacecraft; the x w ~axis is tangent to the 

flight path (relative velocity) and is positive forward; the axes y w ,z w 

are orthogonal to the x w ~axis and are such that the trihedral Px^w^w 1S 

right-handed. In particular, the z w ~axis is contained in the plane of 

symmetry of the spacecraft and is positive downward for the normal flight 

attitude of the spacecraft; the y -axis is positive rightward for the normal 

w 

flight attitude of the spacecraft. The symbols i w »j w » k w denote the unit 
vectors of the wind axes system. 

3 . 5 . Body Axes System . The body axes system Px b y b z b is a Cartesian 
reference frame defined as follows. Its origin P is identical with 

the instantaneous position of the spacecraft; the y b ~axis is orthogonal to 

the plane of symmetry of the spacecraft and is positive rightward; the 

axes x. ,z. are orthogonal to the y, -axis, are contained in the plane of 
b b D 

synwetry, and are such that the trihedral Px b y b z b ^ right-handed. In 
particular, the x b ~axis is positive forward, the y b ~axis is positive 
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rightward, and the z b -axis is positive downward for the normal flight 
attitude of the spacecraft. The symbols ^ b *^b* ^ eno ^ e the unit vectors 
of the body axes system. 



IfHff 

iiijjjiiiii 
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4. Relations between Coordinate Systems 

In this section, the relationships between the different coordinate 
systems are derived; more specifically, attention is focused on the 
following system pairs: inertial axes-Earth axes; Earth axes-local horizon; 
inertial axes-local horizon; local horizon-wind axes; and wind axes-body axes. 

4.1. Basic Relations . In the inertial axes system, a point P can be 
described via its Cartesian coordinates or via its spherical 

coordinates r^e^,^. Here, r i is the radial distance from the center of 
the Earth; e.. is the longitude, positive eastward; and <J>. is the latitude, 
positive northward. The following relations hold betweenCartesian coordinates 
and spherical coordinates: 


x. = r 1 -cose i costly, 

y. j = r.jSine.j cos^, 

z. = r i sin<b i . 


(la) 

(lb) 

(lc) 


In the Earth axes system, a point P can be described via its Cartesian 
coordinates x g ,y e ,z e or via its spherical coordinates r,0,4>. Here, r is 
the radial distance from the center of the Earth; 6 is the longitude, positive 
eastward; and <J> is the latitude, positive northward. The following relations 
hold between Cartesian coordinates and spherical coordinates: 


x g = rcose cos<{>, 


y = rsinQ cos<J> , 


(2a) 

(2b) 


z g = rsin4>. 


(2c) 
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Note that the inertial axes system and the Earth axes system coincide at time 
instant t = 0. Also note that the Earth axes system rotates with constant 
angular velocity w with respect to the inertial axes system. Therefore, 
in spherical coordinates, the following transformation relations hold: 


e i = e + wt, 

4>.j = <f>> 

implying the following inverse relations: 


(3a) 

(3b) 

(3c) 


0 = 0^ “ (JL)t ) 

<J> = 

4.2. Transformation from 
system 0* e y e z e obtained from 


(4a) 

(4b) 

(4c) 


Inertial Axes to Earth Axes . The Earth axes 
the inertial axes system Ox^y^ z - by means 


of the counterclockwise rotation wt around the z^-axis. Note that the 


z -axis is the same as the z.-axis. 
e i 

In vector-matrix notation, the transformation leading from the inertial 


axes to the Earth axes can be expressed as follows: 


cos(wt) 

sin(wt) 

0 


- 

u 






sin(wt) 

cos(wt) 

0 


j i 

-V 

0 

0 

1 


L k i 1 


(5a) 


<3r'i 
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with the implication that 

(5b) 

4.3. Transformation from Earth Axes to Local Horizon . The local 
horizon system Px h y h z h can be obtained from the Earth axes system 0x e y g z e 
by means of the combination of four rotations and one translation. This 
requires the definition of four intermediate coordinate systems, the system 
0Xiy-|Z^; the system 0x 2 y 2 z 2 ; the system Px 3 y 3 z 3 ; an d the system Px^y^z^. 

The system Ox^Zj is obtained from the Earth axes system 0x g y e z e by 
means of the counterclockwise rotation 0 around the z g -axis. Note that the 
z^-axis is the same as the z e ~axis, that the axes x-|,y-| are contained in the 
equatorial plane, and that the axes x 1 ,z ] are contained in a meridian plane. 
The symbols denote the unit vectors of the system Ox^z-j. 

The system 0x 2 y 2 z 2 is obtained from the system Ox^^^ by means of 
the clockwise rotation cj) around the y-j-axis. Note that the y 2 ~axis is the 
same as the y^axis, that the axes x 2 »z 2 are contained in a meridian 
plane, and that the axes y 2 ,z 2 are contained in a plane parallel to the local 
horizon plane. The symbols i 2 »J 2 »k 2 denote the unit vectors of the system 
Ox 2 y 2 z 2 ■ 

The system Px 3 y 3 z 3 is obtained from the system Ox 2 y 2 z 2 by means of 
the radial translation r, leading from point 0 to point P. Since there 
is no rotation, the axes x 3 »y 3 »z 3 are parallel to the axes x 2 ,y 2 ,z 2 » 

are contained in a meridian plane, while 


k i 


cos(wt) -sin(wt) 0 

sin(cot) cos(cot) 0 

0 0 1 


in particular, the axes x 3> z 3 
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the axes y~,z- are contained in the local horizon plane. The symbols 
? 3 , j 3 , k 3 denote the unit vectors of the system Px 3 y 3 z 3 . 

The system Px 4 y 4 z 4 is obtained from the system P* 3 y 3 z 3 by means of 
the counterclockwise rotation ir/2 around the z 3 ~axis. Note that the 
z 4 -axis is the same as the z 3 -axis, that the axes y 4 ,z 4 are contained in a 
meridian plane, while the axes x 4 ,z 4 are contained in the local horizon plane. 
The symbols ^ 4 > J 4 » ^ 4 denote the unit vectors of the system P* 4 y 4 z 4> 

The local horizon system P* h y h z h is obtained from the system Px 4 y 4 z 4 
by means of the clockwise rotation tt/ 2 around the x 4 ~axis. Note that the 
x h -axis is the same as the x 4 -axis, that the axes y h » z h are contained in a 
meridian plane, and the axes x h ,y h are contained in the local horizon plane. 

In vector-matrix notation, the successive transformations leading from 
one coordinate system to another can be expressed as follows: 


■ -+ 

« 

COS0 

sin0 

0 




0 

j l 

- 

-sin0 

COS0 





k l . 


0 

0 

1 










*> “ 

*2 


COS<|> 

0 

sintj) 

-V 




0 

j 2 

— 

0 

1 

1 

ro 


-sin<j> 

0 

COS(|) 



( 6 a) 


( 6 b) 


and 
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Equations (6) imply that 


COS0 coscj) 

sine cost}) 

sin<J> 

-sine 

cose 

0 

-cose sin<j> 

-sin6 sin<J> 

C0S<1> 



while Eqs. (7) imply that 



(7a) 


(7b) 


(7c) 


(8) 


( 9 ) 
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Therefore, upon combining Eqs. (8)-(9), we see that 

(10a) 

with the implication that 

(10b) 

4.4, Transformation from Inertial Axes to Local Horizon . The 
transformation leading from the inertial axes system to the local horizon 
system requires the combination of five rotations and one translation, as one 
surmises from Sections 4.2 and 4.3. If one combines Eqs. (5a) and (10a) 
and accounts for Eqs. (3), the following result is obtained: 


■ ~y 


-sin0 

cos0 si nd> 

-COS0 cos<f> 

Je 

= 

COS0 

sine sin<j> 

-sine cos<}> 

V 

. « 

0 

-cos<J> 

-sin<|> 


■ -y 


-sine 

cose 

0 

^h 

= 

cos0 sintj) 

sine si n4> 

-COS(j> 

SI 


-cose cos<(> 

-sine cos4> 

-sin«J> 


i h 

-y 


- 



■ 


• -y 

-sine.. 


cos0 i 

0 


u 

■y 

cose i 

sin<j> i 

sin6 i si n4> 1 

-C0S(|) j 


j i 

-y 

-COS0. 

cos<}>.. 

-sine,. cos<J>. 

-si ncf> - 


k i 1 


(11a) 


with the implication that 
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(11b) 

4 . 5 . Transformation from Local Horizon to Wind Axes . The wind axes 
system ^ >x w y^ z w can be obtained from the local horizon system Px^y^z^ by 
means of the combination of three rotations. This requires the definition 
of two intermediate coordinate systems: the system PXgy 5 z 5 an< i ^hie system 

Px 6 y 6 z 6* 

The system P^gZg is obtained from the local horizon system P* h y h z h 

by means of the counterclockwise rotation x around the z h ~axis. Note 

that the z^-axis is the same as the z^-axis, that the axes Xg,yg are contained 

in the local horizon plane, and that the axes x 5> z 5 are contained in the 

plane (0P,V), where OP is the radius vector connecting the points 0 and P 

and V is the relative velocity vector, namely, the velocity of the spacecraft 

with respect to the Earth axes system. Also note that the axis x 5 has the 

+ . 

direction of the projected relative velocity vector Vp; this is the projection 
of V on the local horizon. The angle x is called the heading angle and is 
positive if the projected relative velocity vector Vp is directed outward with 

-4 -V -4 

respect to the local parallel. The symbols ^* 35*^5 denote the unit vectors 
of the system PXgyg^* 

The system P x 6 yg z 6 is obtained from the system Px 5 y 5 z 5 by means of 
the counterclockwise rotation y around the y^-axis. Note that the y^-axis 
is the same as the y^-axis and that the axes Xg,Zg are contained in the 
plane (0P,V). Also note that the Xg-axis is positive forward and that the 


4 

u 

j i 

It 

1 

-4 

k i . 



-sin9 i cose i sincfK 

cos0 i s i n 0 ^ sincp^ 

0 -cos<t>.j 


-COS0.j COScj)^ 

-sine.. cos<t>.. 
-sin^ 


-> 

iu 
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Zc-axis is positive downward. The angle y is called the path inclination 
6 -*• 
and is positive if the relative velocity vector V is inclined upward with 

respect to the local horizon. The symbols denote the unit vectors 

of the system Px g y g Zg. 

The wind axes system Px^z w is obtained from the system PxgygZg 

by means of the counterclockwise rotation y around the Xg-axis. Note that 

the x -axis is the same as the x.-axis. Also, note that the x-axis is 
w ° w 

positive forward, the y w -axis is positive rightward, and the z w -axis is 
positive downward and is contained in the plane of symmetry of the spacecraft. 
The angle y is called the angle of bank and is positive if the spacecraft 
is banked to the right. 

In vector-matrix notation, the successive transformations leading 
from one coordinate system to another can be expressed as follows: 


■ 




T 

i 5 


cosy 

siny 

0 

-+■ 


-siny 


0 


— 

cosy 

-* 




1 

. k 5. 


0 

0 






' -> 




* 



cosy 

0 

-siny 







= 

0 

1 

0 

. ^ . 


siny 

0 

cosy 








_ 



i 

w 


1 

0 

0 

j w 

= 

0 

cosy 

siny 




-siny 


1 k w. 

. 

0 

cosy 


k L 


l'r 


(12a) 


(12b) 


(12c) 


I? iff *! 1 
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Equations (12) lead to 


cosy cosy 


cosy smx 


siny siny cosy siny siny siny 

-cosy siny +cosy cosy 

cosy siny cosx cosy siny siny 

+siny siny -siny cosy 


-siny 


siny cosy 


cosy cosy 


with the implication that 


cosy cosy 


cosy siny 


siny siny cosy 
-cosy siny 

siny siny siny 
+cosy cosy 

siny cosy 


cosy siny cosy 
+siny siny 

cosy siny siny 
-siny cosy 

cosy cosy 


4.6. Transformation from Wind Axes to Body Axes . The body system 

Px b y b z b can ol:)ta '' necl ^ rom t * ie W1 ' nc * axes s y stem PVw z w means the 
combination of two rotations. This requires the definition of one intermediate 

coordinate system, the system Px 7 y 7 z 7 . 

The system Px 7 y ? z 7 is obtained from the wind axes system Px^^ by 
means of the counterclockwise rotation o around the z^-axis. Note that the 
z ? -axis is the same as the z w -axis and that the axes x y ,z y are contained 
in the plane of symmetry of the spacecraft. Also note that the axis x 7 
is positive forward, the axis y y is positive rightward, and the axis z 7 
is positive downward. The angle o is called the sideslip angle and is 
positive if the relative velocity vector V is directed leftward with respect 
to the plane of symmetry of the spacecraft. 
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The body axes system Px b y b z b is obtained from the system Px y y 7 z y by 
means of the counterclockwise rotation a around the y 7 ~axis. Note that 
the y b -axis is the same as the y 7 -axis and that the axes x b »z b are contained 
in the plane of symmetry of the spacecraft. Also note that the axis x b 
is positive forward, the axis y b is positive rightward, and the axis z fa 
is positive downward. The angle a is called the angle of attack and is 
positive if the relative velocity vector V is directed downward with respect 

to the x b -axis of the spacecraft. 

In vector-matrix notation, the successive transformations leading 
from one coordinate system to another can be expressed as follows: 


-y 


F 


■ 


• -* 

1 7 


cosa 

sina 

0 


1 w 

J 7 


-sina 

cosa 

0 


j w 

k 7 . 


0 

. 

0 

1 


. V. 


-V " 




■ 


■ -> T 

'b 

-+ 


cosa 

0 

-sina 


h 






-y 



0 

1 

0 


J 7 

-y 







k b J 


sina 

0 

cosa 


. k 7 . 


with the implication that 


■ -V 

'b 

j b 


cosa 

cosa 

cosa sina 

-sina 

= 

-sina 


cosa 

0 

,. k b. 


sina 

cosa 

sina sina 

cosa 

. 

• 






(14a) 


(14b) 


(15a) 
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and that 


■ -y- 


r 



i 


cosa 

cosa 

-sina 

w 





-y 





J w 

- 

cosa 

sina 

cosa 

-> 




0 

k 


-sina 


w 

X • 


. 





(15b) 
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5. Angular Velocity 

In this section, we compute the angular velocity (or evolutory velocity) 
of the local horizon system with respect to the inertial axes system. To 
do so, consider the behavior of the spacecraft between the time instants 
t and t + dt, and denote by d^. the infinitesimal vectorial rotation of 
the local horizon system with respect to the inertial axes system. This 
infinitesimal vectorial rotation can be decomposed into partial rotations 
as follows: 

d “h1 = ®h4 + d “43 + d «32 + d "21 + d "le + di V 

with 

= 

d^43 = 0 , 

dfl 32 = 0, 
d^i — — dcj) j *| , 

dn le = de k e , 

4 ■+ 

dft . = todt k . . 
ei i 

Here, de denotes the infinitesimal change of the longitude, dcf> denotes the 
infinitesimal change of the latitude, and o>dt denotes the infinitesimal 
rotation of the Earth axes system with respect to the inertial axes system. 
The rotation wdt occurs around the z.-axis and is positive counterclockwise; 


(16) 

(17a) 

(17b) 

(17c) 

(17d) 

(17e) 

(17f) 
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the rotation de occurs around the z g -axis and is positive counterclockwise; 
and the rotation d<j> occurs around the y-j-axis and is positive clockwise. 
This explains the difference in the signs appearing on the right-hand 
sides of Eqs. (17d), (17e), (17f). 

Upon combining Eqs. (16)-(17), we see that the infinitesimal vectorial 
rotation of the local horizon system with respect to the inertial axes 
system can be written as 

dft hi = wdt k. + de k g - d<j> j.j. 

As a consequence, the angular velocity of the local horizon system with 
respect to the inertial axes system is given by 

-V ■> —J- 

u hi = dn hi^ dt = uk i + 0k e ‘ 

In the light of Eqs. (5)-(ll), the unit vectors k^ k 0 , j-| can be expressed 
in terms of the unit vectors of the local horizon system as follows: 

k. = -cos<f>. j h - sin<(> i k h , 

-> -y 

k g = -costJ> j h - si nc}) k h , 

- v 

so that 

^hi h - ^ cos ^ + uC0S( t > j ) ” (8sin<j> + usin 4 >^)k^. 


(18) 

(19) 

(20a) 

(20b) 

(20c) 

(21) 


Next, we invoke Eqs. (3), 
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0.j = 0 + ut, 

= 4 >» 

and observe that 


. (22a) 
(22b) 
(22c) 


• • 

0- = 0 + co, 

• * 

4> n - = <f>* 

As a consequence, the angular velocity (21) can be rewritten as 


(23a) 

(23b) 

(23c) 



- e, cos*, j h 


e i 


si ncj>^ 



Next, Poisson's formulas are employed to compute the derivatives of 
the unit vectors of the local horizon system with respect to time: 


(24) 


oydt = : h1 x T h , 

d V dt = “hi x V 

dk h /dt = ; h1 4, 

Upon combining Eqs. (24)-(25), we obtain the relations 
di h /dt =~(e i sincj> i )j h + (e i cosd> i ) k h , 


(25a) 

(25b) 

(25c) 


(26a) 
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dj h /dt = (e i sin^ i )i h - <J> i k h , 
dk h /dt =-(e i cos4> i )i h + 
whose vector-matrix form is the following: 


(26b) 

(26c) 



It is interesting to note that Eq. (27) can also be obtained by taking the 
time derivative of Eq. (11a) and using Eq. (lib). 



23 


AAR-239 


6. Kinematical Equations 

In this section, we derive the scalar relationships corresponding to 
the vectorial equation 

dOP/dt = V.. ( 28 ) 

Here, V. denotes the velocity of the spacecraft with respect to the inertial 
axes system and OP denotes the position vector joining the center of the 
Earth 0 with the spacecraft position P. 

-y 

First, we observe that the position vector OP is given by 

OP = -r.k h , < 29 ) 

-y 

where r. is the radial distance from the center of the Earth and k h is 
the third unit vector of the local horizon system. As a consequence, the 
time derivative of OP can be written as 


dOP/dt = -r.k h - r.(dk h /dt), 

where dk h /dt is given by Eq. (26c). Therefore, upon combining Eqs. (26c) 
and (30), we obtain the relation 


dOP/dt = © i r i cos<f) .i h - <t> i r i - r i k h . 


(31) 


We recall that the relative velocity vector V (velocity of the spacecraft 
with respect to the Earth axes system) can be identified via three elements: 
the relative velocity modulus V, the relative path inclination y, and 
the relative heading angle x ( see Section 4.5). Analogously, the inertial 
velocity vector V i (velocity of the spacecraft with respect to the inertial 
axes system) can be identified via three elements: the inertial velocity 
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modulus V.j, the inertial path inclination and the inertial heading 
angle x r 

Let V p and V ip denote the projections of the vectors V and V i on 
the local horizon. The relative path inclination y is the angle between 
the vectors V and V p ; analogously, the inertial path inclination y i is 

-V 

the angle between the vectors V.. and The relative heading angle x 

is the angle between the vectors V p and i^; analogously, the inertial heading 

•> 

angle x-j is the angle between the vectors V.. p and i^. The conventions used 
for the signs of y i , x-j are analogous to the conventions used for the signs 

of y , x- 

With the above understanding, just as the relative velocity vector V 
can be written in terms of its components on the local horizon system, 

V = Vcosy cosy i h + V cosy sinx j h - Vsiny k h> 

-> 

the inertial velocity vector V i can be written in terms of its components 
on the local horizon system, 

-> -> 

V i = V i cosy i cosx i i h + V i cosy i sinx i j h - V i siny i k h » 

Finally, upon combining Eqs. (23), (31), (33), and upon projecting 
the resulting vectorial equation on the axes of the local horizon system, 
we obtain the following scalar form of the kinematical equations: 

0 . = V 1 -cosy i cosx i /r i cos<j> i , 

J. =-V 1 .cosy i sinx i /r^ , 


(32) 

(33) 


(34a) 

(34b) 


r i = V.siny i . 


(34c) 
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6.1. Relations between Inertial Velocity and Relative Velocity . We 
employ the theorem of composition of velocities, which states that the inertial 
velocity V- is the sum of the relative velocity V and the transport velocity 

v 

V, - V ♦ V t . (35) 


with 



u x OP. 


(36) 


Owing to the fact that 
-> 

w = -wcos<j> j h - usintf) k h = -tocost^ j h - o>sin<t>. k h , 
0 *p - -rk h - -r,k h . 


(37a) 

(37b) 


the transport velocity (36) can be rewritten as 

= u>rcos<j> 1^ = ojr_jCos<t>.j i^. ( 38 ) 

Therefore, upon combining Eqs. (32), (33), (35), (38), and upon projecting 
the resulting vectorial equation on the axes of the local horizon system, 
the following scalar equations are obtained: 

(39a) 
(39b) 


V_jCOsy.j cosx.j - wr.j cos4>^ = Vcosy cosx» 

V.-cosy.j sinx^ - Vcosy sinx> 

VjSiny. = Vsiny. 


(39c) 
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\ i 

b- == 


% j 


Laborious manipulations, omitted for the sake of brevity, lead to the 
following transformation relations: 

? 2 
V. = /[v + 2wrVcosy cosy cos<t> + (wrcos<J>) ] , 

2 2 

tany i = Vsiny//[(Vcosy) + 2ajrVcosy cosy cos4> + (u)rcos<t>) ], 

tany^ = Vcosy siny/(Vcosy cosy + wrcos<}>). 

Equations (40) imply the following inverse relations: 

2 2 
V = /[V^ - 2ojr 1 .V i cosy i cosy i cosef^. + (u)r i cos<t> i ) ], 

2 2 

tany = V i siny i //[(V i cosy. j ) - 2wr i V i cosy i cosy i cos^ + (w^cosc^) ], 

tany = V i cosy i siny i /(V. j cosy i cosy^ - w^cos^.) . 


(40a) 

(40b) 

(40c) 

(41a) 

(41b) 

(41c) 
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7. Dynamical Equations 

In this section, we derive the scalar relationships corresponding 
to the vectorial equation 

T + A + W = ma. = m(dV./dt), ( 42 ) 

where T is the thrust, A is the aerodynamic force, W is the gravitational 
force, m is the mass of the spacecraft, and a. is the inertial acceleration. 

We consider the case where the engine is shut-off, so that 


and the mass of the spacecraft is constant. Hence, Eq. (42) is written as 

A + W = ma.j = m(dV^/dt). ( 44 ) 

We now compute the components of the vectors appearing in Eq. (44) on the axes of 
the local horizon system. 

7,1, Aerodynamic Force . The components of the aerodynamic force on 
the wind axes are the drag D, the side force Q, and the lift L. Therefore, 
the aerodynamic force can be written as 



No special significance is implied in the signs appearing on the right-hand 
side of Eq. (45). These signs merely reflect the conventions adopted in 
this report with regard to the positive values for the drag, the side force. 


( 45 ) 


and the lift. 
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In vector-matrix form, Eq. (45) can be rewritten as follows: 


D 0 0 
0 Q 0 
0 0 L 



(46) 


Therefore, upon combining Eqs. (13a) and (46), we obtain the relation 



-Dcosy cosx -Dcosy sinx Dsiny 


’ T 


-Qsiny siny cosy -Qsiny siny siny -Qsiny cosy 



A = 

+Qcosy siny -Qcosy cosy 

i 




-Lcosy siny cosy -Lcosy siny siny -Lcosy cosy 


^h 


-Lsiny siny +Lsiny cosy 




(47) 


The next step consists of rewriting the elements of the matrix in 
Eq. (47) in terms of inertial velocity elements V. ,y. instead of 
relative velocity elements V,y,y • For this purpose, let the following 
functions of the inertial velocity elements be defined: 


M-| = M-j (V.j ,y.j ,x.j ) = V i cosy i cosx^ - wr-cos^, 

M 2 = M 2 ( v i »Yi » x i) = V.cosy i sinx i , 

M 3 = M 3 (V 1 -,Y i .X,) ■ v^iny,, 

N, = NjfVpy^Xj) ■ /(Mi + m| + M3) , 

N 2 = NjfV^y^x,) = /(«i + M|)//(Mf + * m|). 


(48a) 

(48b) 

(48c) 


(49a) 


(49b) 
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N 3 = N 3 (V i ,Y.,x i ) = M 3 //( M 2 + + M 3 ), (49c) 

N 4 = N 4 (V.,Y i »X i ) = + M^), (49d) 

N 5 = N^V-.Y^X-j) = M 2 //( M 2 + M^). (49e) 

In the light of Eqs. (39), the following relations can be established: 

M-j = Vcosy cosx, (50a) 

M 2 = Vcosy sinx, (50b) 

M 3 = Vsinx, (50c) 

and 

N 1 = V, (51a) 

N 2 = cosy. (51b) 

N 3 = sinY, (51c) 

N 4 = cosx, (51d) 

N 5 = sinx. (51e) 


Therefore, upon combining Eqs. (47) and (50)-(51), we obtain the relation 


-DN 2 N 4 

-DN 2 N 5 

dn 3 


’ 4 

1 h 

-Qsiny N 3 N 4 

-Qsiny N 3 Ng 

-Qsiny N 2 


^h 

+Qcosy Ng 

-Qcosy N 4 




-Lcosy N 3 N 4 

-Lcosy N 3 Ng 

-Lcosy N 0 

C 


*h 

-Lsiny N g 

+Lsiny N 4 



. . 
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7.2. Gravitational Force . Here, we assume that the Earth is an 
oblate planet and that its mass has radial symmetry with respect to the 
axis of rotation. Because the equatorial radius r g is larger than the polar 
radius r , the gravity force W has two components: the radial component mg, 
directed toward the center of the Earth, and the latitudinal component mf, 
tangent to the local meridian and directed toward the Equator. Therefore, 
the gravity force can be written as 

W = mfj h + mgk h . 

The radial component g and the latitudinal component f of the 
acceleration of gravity are related to the Earth's gravitational potential 
U by the expressions 


g = 3 U/ 3 r.. , f = (1/r.. )3U/3<J>.. , 


(54) 


where 


U Myrpn + J^r-g/r,) 2 ^ + J 3 (r e /r- i ) 3 H 3 + > 

H 2 = 1/2 - (3/2)sin 2 <J>., 

H 3 = ( 3 / 2 )sin<f>. - ( 5 / 2 )sin 3 <f>., 

H 4 =-(3/8) + (30/8)sin 2 4>. - (35/8) si n 4 ^. . 


(55a) 

(55b) 

(55c) 

(55d) 


Here,y g is the Earth's gravitational constant, r g is the equatorial radius, 
and J 2 » J 3 » J 4 denote the characteristic constants of the Earth's 
gravitational field. Note that the expression for U is approximate, since 
harmonics of order higher than four are ignored. 
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Upon combining Eqs. (54)-(55), we see that the components of the 
acceleration of gravity can be written as 

g * (u e /rf)[l +J 2 ( r e /r i )2G 2 + J 3 (r e /r i )3G 3 + J 4 (r e /r i 'V ’ 

G 2 = 3/2 - (9/2)sin 2 (j) i , 

3 

G 3 = 6si n^^ - lOsin , 

G 4 =-l 5/8 + ( 150/8) si n 2 <f> i - ( 175/8) si n 4 ^ , 


(56a) 

(56b) 

(56c) 

(56d) 


and 


f ■(u e /r^)[J 2 (r e / r 1 ) 2 F 2 + J 3 (r e /r,) 3 F 3 + )\l , 

F 2 = 3sin4>.. cos<J>.j , 

F 3 =-(3/2)cos(j>. + (15/2)sin 2 <J> i cos4> Jj , 

F 4 =- ( 1 5/2 ) si nd> i cos<J>. + ( 35/2 ) si n 3 cj> i cos^. 


(57a) 

(57b) 

(57c) 

(57d) 


7.3. Inertial Acceleration . Let V ixh , V. yh , V- zh denote the components 
of the inertial velocity on the local horizon system. 


V 


ixh 



= V.cosy-j cosx i5 
= V i cosy i sinx i > 
=-V.sin Yi . 


With this understanding, the inertial velocity (33) can be rewritten as 


(58a) 

(58b) 

(58c) 


V i ” V ixh n h + V iyh J h + V izh k h* 


(59) 
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Therefore, the inertial acceleration is given by 


d V dt ■ v 1xh 1 h + V iyh j h + v 1zh k h 

+ V ixh (d V dt) + v iyh (d V dt) + V izh (d V dt) ‘ 

If we combine Eqs. (26) and (34), the time derivatives of the unit vectors 
of the local horizon system can be written as 

d? h /dt = -(V i cosy i cos Xi tan<j>./r.)j h + (V.cosy i cos X ./r.)k h , 

dj h /dt = (V^cosy.j cos Xi tan<J> i /r i )i h + (V i cosy i sin^./r.. )k h , 

dk h /dt = -(V i cosy i cos Xi /r i )i h - (V i cosy i siny./^)^. 

Upon combining Eqs. (58), (60), (61), the inertial acceleration becomes 

dV^/dt = (V ixh + V^cos 2 y i cos Xi siny. tarw^/r. + V?cosy i sin Xi cos Xi /r.)i h 

+ (V iyh - V 2 cos 2 y i cos 2 y. tan^/r. + V 2 cosy i sin Xi sin x ./r i )j h 

+ < J izh + vfcos 2 y 1 /r 1 )k h . 


(60) 


(61a) 

(61b) 

(61c) 


(62) 


7.4. Scalar Equations . Next, we combine Eqs. (44), (52), (53), (62). 
Upon projecting the resulting vectorial equation on the axes of the local 
horizon system, we obtain the following scalar equations: 


V ixh = -DN 2 N 4 /m + (QN 5 /m)cosy - (QN 3 N 4 /m)siny 

- (LNg/m)siny - (LN 3 N 4 /m)cosy 

- (V 2 /r i )(cos 2 y i cos Xi sin Xi tan<t> i + cosy i siny i cosy..). 


(63a) 
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V. yh = -DN 2 N 5 /m - (QN 4 /m)cosy - (QN 3 N 5 /m)siny 
+ (LN 4 /m)siny - (LN 3 N^/m)cosy 

0 o O 

+ f + (V^/r^Ccos Y i cos Xi tan<|>. - cosy i siny i sinx.,0, (63b) 

V izh = DN 3 /m - (QN 2 /m)siny - (LN 2 /m)cosy 

+ g - ( V?/r^ )cos^y 1 - . (63c) 

We recall that the components of the inertial velocity on the axes of the 
local horizon system are given by [see Eqs. (58)] 


V ixh = V i C0SY i cosx i’ 
V iyh = V i C0SY i sinx i * 
V izh = - V i sin V 


with the implication that 


(64a) 

(64b) 

(64c) 


>h 


cosy.. cosy^ 

-siny i cosx i 

-sinx i 


• 

V i 


V iyh 

= 

cosy. sinx i 

*siny. sinx i 

cosx i 


Vi 

, (65) 



-siny i 

-cosy.. 

0 


V i cosy i x i 



and that 


T ♦ " 

v i 


cosy i 

cos Xi 

cosy^ 

sinx i 

-siny. 


1 

sz 

X 

• > • 

■ 

Vi 

= 

-siny i 

cos Xi 

-siny i 

sinx i 

-cosy. 


V iyh 

• 

• 

V i cosy i x i 


-sinx i 


cosx i 


0 


. v i* .. 


( 66 ) 
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The final step consists of combining Eqs. (63) and (66). This leads to 
the following scalar form of the dynamical equations: 


V. = -DN 7 /m + (DN 6 /m)cosy i cos X .. 


- (QN 4 N 6 /m)siny siny. + (QNg/mN 2 )cosu cosy i sin Xi 


- (LN 4 N 6 /m)cosy siny i - (LNg/mN 2 )siny cosy i sin Xi 


+ fcosy^. sin Xi - gsiny^ , 


(67a) 


V i Y i =-(DN 6 /m)siny i cos Xi 

+ (Q/mN 2 )siny cos Yi - (QNgN 7 /mN 2 )siny sin 2 Yi cos X . - (QNg/mN 2 )cosy sin Yi sin Xi 
+ (L/mN 2 )cosy cosy i - (LNgN 7 /mN 2 )cosy sin 2 y i cos X . + (LNg/mN 2 )siny siny i sin Xi 


- fsiny i sin Xi + (V 2 /r i - g)cosy i5 


(67b) 


V i cosY iXi =-(DNg/m)sin Xi 

- (QNgN 7 /mN 2 )siny siny^ sin X ^ - (QN 7 /mN 2 )cosy cosy.. + (QNg/mN 2 )cosy cos Xl - 

- (LNgN 7 /mN 2 )cosy siny. sin Xi + (LN 7 /mN 2 )siny cosy i - (LNg/mN 2 )siny cos Xi 

(67c) 


? 2 

+ fcos Xi + (V./r^cos y i cos Xi tan#. 


where 


ig = wr i cos(|)./N 1 , 


N 7 = V N 1 


(68a) 

(68b) 


Equations (67) can be rewritten as 
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V.j =-DN 7 /m + (DNg/m)cosY 1 . cosx^ 

- (QN^Ng/m)siny siny. + (QNg/mN 2 )cosy cosy i sinx i 

- (LN^Ng/m)cosvi siny i - (LNg/mN 2 )siny cosy.. sinx i 

+ fcosy i sinx i - gsiny^, ( 69a ) 

Yi =-(DNg/mV i )siny i cosx i 

+ (Q/tnN 2 V.)siny cos Yi - (QNgN 7 /mN 2 V . )siny sin 2 Yi cos Xi - (QNg/mN 2 V i )cosy siny. sin X . 

+ (L/mN 2 V.)cosy cos Yi - (LNgN 7 /mN 2 V. )cosy sin 2 Yi cos Xi + (LNg/mN^ )s1ny siny. sin X . 

- (f/V i )siny i siny^ + ( v ^/ r ^ “ g/V^Jcosy.., (69b) 

X-j =-(DNg/mV i )sinx i /cosy i 

- (QNgN 7 /mN 2 V i )siny tan Yi sin Xi - (QiymNgV. )cosy + (QNg/mNgV. )cosy cosx 1 -/cosy i 

- (LNgN 7 /mN 2 V 1 - )cosy tan Yi sinx^ + (LN 7 /mN 2 V. )siny - (LNg/mN 2 V^ )siny cosx^/cosy^ 

+ (f/V.)cosx i /cosy i + (V i /r i )cosy i cosx i tan<J>.. 


(69c) 
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8. Summary of Results 

In this report, we have derived the equations of motion of a spacecraft 
under the following assumptions: (a) the spacecraft is a particle of constant 
mass; (b) the Earth is rotating with constant angular velocity; (c) the 
atmosphere is at rest with respect to the Earth; (d) the Earth is an oblate 
planet, and the gravitational potential depends on both the radial distance 
and the latitude; however, harmonics of order higher than four are ignored. 

An inertial axes system has been used, and the following kinematical 
and dynamical equations have been obtained: 


• 

6^ = V i cosy i cosx^/r^cosef)^ , 

(70a) 

l i =-V i cosy i sin Xi /r i , 

(70b) 

r i = V^siny. , 

(70c) 

-DN 7 /m + (DN g /m)cosy 1 . cosx i 


(QN^N g /m)siny siny.. + (QNg/mN 2 )cosy cosy^ sinx. 


(LN 4 Ng/m)cosy siny. - (LNg/mN 2 )siny cosy i sinx i 


fcosy^ sinx i - gsiny^ , 

(71a) 


y. =-(DN 6 /mV.)sinY. cosx i 

+ (Q/mty/^sinp cosy i - (QNg^/mN^ )siny sin 2 y i cosx i - (QNg/mfi^ )cosy siny. siny 
+ (L/mN^Ocosy cosy. - (ENgf^/mN^ )cosy sin^ cosx i + (l-Ng/mN^ )siny siny i sinx 
- (f/V i )siny i sinx i + (V i /r i - g/V^cosy^ 


(71b) 
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Xi =-(DN 6 /mV i )sinx i /cosy i 

- (QNg N 7 /mN 2 V.)siny tany i sin Xi - (Q^/mN^ )cosy + (QN 6 /mN 2 V i )cosy cosx i /cosy i 

- (LN 6 N 7 /mN 2 V i )cosy tany i sinx i + (Lf^/mN^ )siny - (LN 6 /mN 2 V i )siny cosx 1 /cosy i 

+ (f/V i )cosx i /cosy i + (V i /r i )cosy i cosx i tan^. (71c) 

In the dynamical equations (71), the quantities N.. depend on the inertial 
velocity elements , y^ , x-j fl nd are given by 


N-j = /(M* + M 2 + M 3 ), 

(72a) 

N 2 = /(M^ + M 2 )//(M^ + M 2 + M 3 ), 

(72b) 

n 3 = m 3 //(m^ + m 2 + M 3 ), 

(72c) 

N 4 = M-j//(M^ + M 2 ) , 

(72d) 

n 5 = m 2 //(m^ + m 2 ). 

(72e) 

Ng = ur^cos^^/N-i , 

(72f) 

N y = V i /N 1 , 

(72g) 

M-j = V i cosy i cosx^ - wr^cosfj)^ , 

(73a) 

M 2 = V.cosy i sin Xi , 

(73b) 

M 3 = V^lny,. 

(73c) 
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8.1. Aerodynamic Force . In Eqs. (71), the drag, the side force, 
and the lift are given by 

n 2 2 2 

D = (1/2)C d pS(V^ - 2wr i V i cosy i cosy^ cos^. + w r^cos 4^), (74a) 

_ 2 2 2 

Q = (1 /2)CgpS(V? - 2wr i V i cosy i cosy i cos<j> i + u> rfcos d> i ) , (74b) 

n 2 2 2 

L = (1/2)C l pS(V^ - 2wr i V i C0SY i cosy i cos<t> i + orr.cos <4> i ) » (74c) 

where C Q is the drag coefficient, Cg is the side force coefficient, C L 
is the lift coefficient, p is the air density, and S is a reference surface 
area. In turn, the aerodynamic coefficients are functions of the form 


C D = C D (a,o,M,R e ) , 
Cq = Cg(a,a ,M,R g ) , 
C L = Cg(a,a,M,R e ) , 


(75a) 

(75b) 

(75c) 


where a is the angle of attack, a is the sideslip angle, M is the Mach 
number, and R g is the Reynolds number. 

8.2. Gravitational Force . In Eqs. (71), the radial component and the 
latitudinal component of the acceleration of gravity are given by 


g - (p e /r?)n + J 2 (r e /i' j ) 2 G 2 + + J 4 (r e / r 1 )\] ■ 

G 2 = 3/2 - (9/2)sin 2 J r 

3 

G 3 = 6sin4> i - lOsin <J> i , 

G 4 =-l 5/8 + (150/8)sin 2 <j). - ( 1 75/8) sin 4 ^^. , 


(76a) 

(76b) 

(76c) 

(76d) 
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and 


f =(y e /r^)[J 2 (r 

e /r,) 2 F 2 + 0 3 (i- e /r 1 ) 3 F 3 + J,(r e /r, )%] . 

(77a) 

F 2 = 3sin<J)^ cos<J) 

1* 

(77b) 

F 3 s= -(3/2)cos4>- 

+ (15/2)sin 2 <J>.. cos<j>.. 

(77c) 

F 4 =-(15/2)sin<j>. 

3 

cos<J> i + ( 35/2 ) sin ^ cos<t>... 

(77d) 

8.3. Physical Constants. The major physical constants appearing in 
the system (70)-(77) have the following values: 


u = 0.729211595 

E-04 rad/sec. 

(78a) 

y e = 0.39860064 

E+15 m 3 /sec 2 , 

(78b) 

J 2 = 0.10826271 

E-02, 

(78c) 

J 3 =-0.25358868 

E-05, 

(78d) 

J 4 =-0.1624618 

E-05, 

(78e) 

r = 0.6378164 
e 

E+07 m. 

(78f ) 

r = 0.6356755 
P 

E+07 m. 

(78g) 

Here, w is the Earth's angular velocity;y e is the Earth's gravitational 
constant; J 2 , J 3 , J 4 are the characteristic constants of the Earth's gravitational 
field; r 0 is the Earth's equatorial radius; and r p is the Earth's polar radius. 
Note that the Earth's sea-level radius r $£ varies with the latitude <|>. 


according to the relation 
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r s£ = (l/2)(r e + r p ) + (l/2)(r e - r p )cos(2(f) i ) . (79) 

8.4. Spacecraft Data . For the AFE vehicle, it is assumed that 
m = 0.16782918 E+04 kg, (80a) 

S = 0.14314 E+02 m 2 , (80b) 

a = 0.17000 E+02 deg, (80c) 

C L =-0.370696 E+00, (80d) 

C D = 0.131452 E+01 . (80e) 


Here, m is the spacecraft mass at atmospheric entry; S is the reference 
surface area; a is the angle of attack; is the lift coefficient; and Cq 
is the drag coefficient. Note that, for the aeroassisted flight experiment, 
the angle of attack is kept constant; the aerodynamic coefficients are assumed 
to be independent of the Mach number and the Reynolds number; and the spacecraft 
is controlled via the angle of bank. 



41 


AAR-239 


9. Transformation Relations 

In this section, we summarize the transformation relations which allow 
one to pass from (i) quantities computed in an Earth-fixed system to 
(ii) quantities computed in an inertial system, and viceversa. 

9.1. Spacecraft Position . Let r, 6, <J> denote the spherical coordinates 
of the spacecraft P in the Earth-fixed system Ox e y e z e . Let r^ , 0^, 
denote the spherical coordinates of the same spacecraft in the inertial 
system Ox^^. Assume that the axes of the Earth-fixed system coincide 
with the axes of the inertial system at time instant t = 0. Then, the 
following transformation relations hold: 


0. = 0 + U)t, 

4>-j = <f>. 

Equations (81) imply the following inverse relations: 


(81a) 

(81b) 

(81c) 


(82a) 


0 - 0. - 0)t, 

4> = 4>-j. 

9.2. Spacecraft Velocity . Let V, y, x denote the velocity modulus, the 
path inclination, and the heading angle in the Earth-fixed system Ox e y g z e . 

Let V. , y. , x-j denote the velocity modulus, the path inclination, and the 
heading angle in the inertial system Ox^y^z.. The following transformation 


(82b) 

(82c) 


relations hold: 
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? 2 

V. = /[V + 2wrVcosy cosy cos<j> + (wr cos<}>) ] , 

? 2 

tany i = Vsiny//[(Vcosy) + 2iorVcosy cosy cos<}> + (u>r cos<J>) ], 

tany^ = Vcosy sinx/(Vcosy cosx + wrcost])). 

Equations (83) imply the following inverse relations: 

V = Av\ - 2cor^V^cosy 1 . cosy^ cos^^ + (wr^coscj)^ )^] , 

tany = V.siny.//[(V.cosy i ) 2 - 2u>r i V. j cosy i cosx i cos(J) i + (accost}). ) ]> 

tanx = V.cosy i sinx i /(V i cosy i cosx i - w^cost}).. ). 


(83a) 

(83b) 

(83c) 

(84a) 

(84b) 

(84c) 


9.3. Cartesian Coordinates . After the spacecraft position is known 
in spherical coordinates, the corresponding Cartesian coordinates can be 
computed. The following transformation relations hold: 


x g = rcose cosef), 
y g = rsinG cost}), 
z g = rsintj), 


(85a) 

(85b) 

(85c) 


and 


= r..cos0.j cos<j) 1 . , 


y.. = r^sine^ cosef)^ , 


z i = r i sinc|> i . 


(86a) 

(86b) 

(86c) 
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9.4. Cartesian Velocity Components . After the spacecraft velocity 
elements V, y, x V.» y.. , y.j are known, the Cartesian velocity components 
x g , y g , z e or , y.j , z. can be computed. The following transformation 
relations hold : 

x g =-Vsin0 cosy cosy + Vcos9 sin<j> cosy siny + VcosQ cos<J> siny, 

y g = V cose cosy cosx + Vsine si n<f. cosy sinx + Vsine cos<}> siny, 

z =-Vcos<f> cosy siny + Vsin<f> siny. 


and 


=-V.sine. 

cosy^ 

cosy i 

+ v i cose i 

cos4>^ 

siny i , 

= v i cose i 

cosy^ 

cos Xi 

+ v^sine^ 

COScp^ 

siny i , 


cosy i 

sin Xi 


(87a) 

(87b) 

(87c) 

(88a) 

(88b) 


(88c) 
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10. Conclusions 

This report is the second of a series dealing with the determination 
of optimal trajectories for the aeroassisted flight experiment (AFE). The 
AFE refers to the study of the free flight of an autonomous spacecraft, 
shuttle-launched and shuttle-recovered. Its purpose is to gather atmospheric 
entry environmental data for use in designing aeroassisted orbital transfer 
vehicles (AOTV). 

It is assumed that: the spacecraft is a particle of constant mass; 
the Earth is rotating with constant angular velocity; the Earth is an 
oblate planet, and the gravitational potential depends on both the radial 
distance and the latitude; however, harmonics of order higher than four are 
ignored; the atmosphere is at rest with respect to the Earth. 

Under the above assumptions, the equations of motion for hypervelocity 
atmospheric flight (which can be used not only for AFE problems, but also for 
AOT problems and space shuttle problems] are derived in an inertial system. 
Transformation relations are supplied which allow one to pass from quantities 
computed in an inertial system to quantities computed in an Earth-fixed system, 
and viceversa. 
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